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The parametrically driven, damped Gross-Pitaevskii equation, which models Bose-Einstein con- 
densates in which the interatomic s— wave scattering length is modulated in time, is shown to support 
spatially modulated states in the form of rolls. A Landau equation with broken phase symmetry is 
shown to govern the dynamics of the roll amplitude. 

I. MODEL 

We consider spontaneous pattern formation in dilute Bose-Einstein condensates (BEC) whose interatomic s— wave 
scattering length is periodically varied in time 0. This modulation can be achieved by different means like, e.g., the 
use of magnetic 0, electric or light B fields. When damping is considered || such BEC can be described by the 
following driven, damped Gross-Pitaevskii (GP) equation, 

dtif> (r, t) = (1 - i 7 ) [-V 2 ^ + V (r) ^ + V] + IVf i>- (1) 

which has been written in appropriate normalized variables. Parameter 7 accounts for damping (values on the 
order 0.01 — 0.1 seem to be appropriate for actual BECs ||), a(t) represents the instantaneous value of the 
interatomic s— wave scattering length and do represents its mean value. We shall assume the simplest case 
a (t) = ao [1 + 2a cos (2ut)] . 

In the following we shall consider "pancake" (2D) or "cigar" (ID) shaped BECs in which the trapping potential V (r) 
strongly confines the condensate in one direction or two directions, respectively, whilst along the other direction(s) it 
extends sufficiently as compared with the typical wavelength of the emerging pattern. If a parabolic trapping potential 
V (r) = —I + 5 (w 2 a; 2 + ujyi/ 2 + cj 2 z 2 ) is used (the arbitrary offset is set to —1 for mathematical convenience) the 
above situation can be fulfilled whenever lo z 3> lu 3> u> x 1 f° r the 2D case and lo v ,uj z 3> w 3> lo x for the ID case. In 
a first approximation these inequalities allow: (i) to neglect the confined direction(s) in the descripition of the BEC 
dynamics, and (ii) to approximate the potential by a constant. Both approximations come from the fact that the 
characteristic wavelength of the selected pattern turns out to be, at the same time, much smaller than the size of the 
weakly confined direction and much larger than the strongly confined direction. Hence we consider along the rest of 
this work the following GP equation (l],|(| : 



d t ijj(x,t) = -d'^- tp + |V>| 2 ^ +2acos(2cjt)\^tp. (2) 

Eq. (0) admits the following homogeneous state 

ip = exp \—i (a/uj) sin (2u)t)] , (3) 

which, in the absence of modulation (a = 0) reduces to the BEC ground state ip = 1 (the chemical potential is null 
in this representation because of the choice of the offset in the trapping potential) . 

II. LINEAR STABILITY ANALYSIS: THE PARAMETRIC RESONANCE 

We wish to know whether the spatially homogeneous external driving is able to induce a spontaneous spatial- 
symmetry breaking of (^). Fot that we perform next a linear stability analysis of @ by adding a small perturbation 
to that solution in the form tf> — exp [—i {a/uj) sin (2u>f)] [1 + w (t) cos (kx)]. Substitution of this expression into Eq. 
(||) and linearization with respect to w leads to the following coupled equations for u = Rew, and v = Imw: 
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u' (t) = -7 (2 + fc 2 ) u + k 2 v, (4) 
v' (t) = -jk 2 v - [2 + k 2 + 4a cos (2ut)] u, (5) 

which can be combined to yield |jj 

u" (t) + 2 7 (l + fc 2 ) u' (t) + [(1 + 7 2 ) « 2 (fc) + 4fc 2 a cos (2a;*)] u (t) = 0, (6) 

where 



fi(Jfe) = fcv^ + fc 2 , (7) 

is the nonlinear dispersion relation for the perturbations. 

Eq. (^) is a Mathieu equation with damping, analogous to that describing parametrically driven, damped pendula, 
and ubiquitous in the description of parametric forcing 0. Its solutions are well known which, according to Floquet's 
Theorem, can be written as 

u(t)=Ref(t)e^, (8) 

where / is a periodic complex function of period ir/u> and /i is the so called Floquet exponent (in the case of the 
Mathieu equation fj,/i is known as Mathieu characteristic exponent). The BEC ground state (||) will be unstable 
against perturbations with wavenumber fc whenever Re/x > 0. A general property of Eq. is that Re/i > within 
a series of resonance " tongues" very much like the parametric resonances observed in liquids vibrated vertically . 
For small 7 and a (the cases we consider here) these tongues are located, as a function of f2, around f2 (fc„) = noj, 
n = 1,2,3,... Within these tongues the BEC ground state undergoes a spontaneous spatial-symmetry breaking 
and perturbations with wavenumber fc„ = y/y/l + n 2 ui 2 — 1 amplify. The Floquet exponents can be numerically 
determined by using standard mathematical methods Useful analytical information can be obtained in the limit 
of weak damping 7 <§; 1. Taking into account that for 7 = a = the solution to (||) is of the form u (t) — Re Ue int , 
if we consider a ~ \u — Q\ ~ 7, a perturbative expression for u can be obtained by allowing U to be a slowly varying 
function of time with the result 

u (t) = Re [U+ exp ( A+f ) + U- exp (A+t)] e int , (9) 
A± = i (o> - SI) - 7 (1 + A: 2 ) ± ^ 2 - (a, - «) 2 . (10) 
Note that (||) can be written as (JsJ) with / (t) = exp (i2o>i), and 

n = -™ - 1 (1 + k 2 ) + ^ (^p) "-(cu-n) 2 . (11) 

Finally the condition Re fj, > reads 

a>^^ 2 (l + k 2 ) 2 + (cj-n)\ (12) 

where the equality defines the boundary (neutral stability line) of the first resonance tongue. Note that ( |l2| ) indicates 
that, for fixed to, the threshold for pattern formation is minimum at O = u>, i.e., for wavenumbers fc = fci — 
VvT+ — 1. The minimum parametric driving is hence predicted to be ajjjfn ( w ) = 7^Vl + a; 2 / (Vl + a; 2 — l) . 
Following a similar analysis for the second resonance tongue (Q ~ 2w) in the limit a ~ \u> — f2| ~ 7 1 ^ 2 , it follows that 
w (t) can be written as (^) with 

1 / /(T 2 fc 4 \ 2 rv 2 fc 4 

A ± =z(2o;~r!)-7(l + fc 2 )±-W20f^ r J - 24-^- (2w - fi) - 9 (2c^ - ft) 2 , (13) 
from which the second tongue (Re A+ > 0) runs: 

fc 2 V 10 



a> '-^rrxli^xl \l^l 2 {^ + k 2 ) 2 + 9(2cj-0) 2 + 2(2w-fi). (14) 



The minimum driving amplitude needed in this case, a^Tn > occurs at ft = 2a;. In this case a™-~ 2 ~ \/7 hence the 
second resonance tongue is excited at larger driving amplitudes than the first one, for which aj^ 1 ~ 7. Note in both 
cases that for 7 = (no damping) the threshold for both (and in fact any) resonance tongues is a — hence all are 
simultaneously excited at vanishingly small values of driving amplitude. Damping hence breaks this degeneracy and 
selects the first resonance tongue at low drivings. 
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III. THE ROLL PATTERN: AMPLITUDE EQUATION 



In the following we study the simplest pattern supported by the GP equation (|2|) under the previously described 
parametric instability. We consider a roll pattern in the form 

ip (x, t) = e~ l{a/u>) sin(2wt) [l + ewi (i) cos (kx) + e 2 w 2 (x, t) + e 3 w 3 {x, t) + O (e 4 )] , (15) 

where < e <C 1 is an auxiliary small parameter. In order to deal with a small roll component (of order e) we assume 
that the amplitude of the parametric driving is also small, say a = O (e 2 ) (that this is the proper scaling for a is 
justified a posteriori by the consistency of the final result). On the other hand for the parametric excitation to be 
effective Eq. (|lj) must be fulfilled. If we want to take into account both the resonance condition and the effect of 
damping we must impose 7 = (e 2 j and (u> — Q) = O (e 2 )- Summarizing we consider in the following the scalings 

u = £1 + e 2 u>2, a = e 2 a 2 , 7 = e 2 72- (16) 

Our goal is to find an equation for the roll complex amplitude w\. This will be done using a standard multiple 
timescale technique |^,^| . For this purpose we introduce a slow time 

t = e 2 t, (17) 

and allow all coefficients of the expansion to depend formally both on t and t: 

wi (t) = un (t, t) + iv n (t, t) , (18) 
Wj (x, t) — Uj (x, t, t) + ivj (x, t, t) . (19) 

Finally Eqs. |[|), (|l|), (|l|) and are introduced into the GP equation (g). After using the chain rule for 
differentiation <9 t — > dt +e 2 d T and equating equal powers in e an infinite hierarchy of differential equations is obtained. 



A. Order e 1 



This is the first nontrivial order and reads 

d t U 11 (t,T)-k 2 V 11 (t 7 T) -0, (20) 

d t vn (t, r) + (2 + k 2 ) u n (t, r) = 0, (21) 

whose solution can be written as 

Ull (t, T ) = [r( T )e int +r( T )e~ mt ] , (22) 

vn(t,r) =l ^[ r ( T )e mt ~r(T)e- mt ], (23) 

where r (r) stands for a yet arbitrary complex function of the slow time, and the overbar denotes complex conjugation. 



B. Order e 2 

At this order we find 



d t u 2 (x, t, t) + d 2 v 2 (x, t, t) — [1 + cos (2kx)] uu (t, r) vu (t, t) , (24) 
d t v 2 {x, t, r) + (2 - d 2 x ) u 2 (x, t,r) = -| [1 + cos (2kx)} [Zu 2 ^ (t, r) + v 2 xl (t, t)] , (25) 



whose solution can be written as 



u 2 (x,t,r) — u 20 (t, t) + u 22 (t, t) cos (2kx) , (26) 
v 2 (x, t, t) = v 20 (t, t) + v 22 (t, t) cos (2kx) , (27) 



where each of the coefficients verifies 
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d t u 20 (t, r) = [r 2 (r) e i20t - r 2 (r) e -' i20t ] , (28) 

ft«2o(*,r) =-2u 2Q (t,r)-2(2 + k- 2 )\r(T)\ 2 + (k- 2 -l) [r 2 (r) e ant r 2 (r) e~ ant ] , (29) 

ft U22 (*, r) = 4fc 2 « 22 (t, r) + [r 2 (r) e* 2nt - r 2 (r) e -' i2nt ] , (30) 

d t v 2 2 (t, r) - - (2 + 4fc 2 ) u 22 (t, r) - 2 (2 + fc~ 2 ) |r (r)| 2 + (fc~ 2 - l) [r 2 (r) e l20t r 2 (r) e - l2nt ] , (31) 
which solved yield 

u 20 (t, r) = -2 (1 + fc- 2 ) |r (r) | 2 + -L [ r ( T ) e «* + f ( T ) e -«*] 2 , (32) 

u 22 (t,r)=~[r (r) e' f2t + f (r) e~ int ] 2 , (33) 

« 20 ft r) = v 200 (r) + ± [r 2 (r) - r 2 (r) e"«*] , (34) 

« 22 (t, r) = - S [r 2 (r) e *» - r 2 (r) e"«*] , (35) 



where W200 ( r ) is a yet undetermined function of the slow time, which is not fixed by the present analysis. This 
information should be obtainable by extending the calculation up to higher orders of the expansion. Anyway, as will 
be seen below, the knowledge of i^oo ( T ) is not relevant for our purposes. 



C. Order e s . The amplitude equation 

This is the last order we will consider. It reads 

d t u 3 (x, t, r) + d 2 x v 3 (x, t, t) = f u (t, r) cos (kx) + g u (t, r) cos (3kx) , 
d t v 3 (x, t, t) + (2 - d 2 ) u 3 (x, t, t) = f v (t, t) cos (kx) + g v (t, t) cos (Zkx) , 



where 



72 (2 + fc 2 ) Mu + (2v 2 o + w 22 ) Mil + (2u 2 o + u 22 ) uii + f ("n + u?i) "u, 

2„ 



fu = -d T un 

f v = —d T vn — "f 2 k 2 vn — 4a 2 cos (2Qt + 2uu 2 t) « u 
9u = unv 22 + u 22 vn + \ (u 2 n + vf-^) v n , 
g v = -3wnu 22 - vnv 22 - \ (u\ x + uj x ) uii. 

The solution to Eqs. pq) and (p7|) can be written as 



3 (2u 20 + m 22 ) "li - (2« 20 + « 22 ) «n - | (uli + v\i) 



u 3 (x, t, t) = u 3 i (t, t) cos (fcc) + U33 (t, t) cos (3kx) . 
«3 (a;, £, t) = U31 (i, t) cos (kx) + v 33 (t, t) cos (Zkx) , 



where each of the coefficients verifies 



d t u 3 i (t, t) - fc 2 w 3 i (t, t) = f u (t, t) , 

d t v 31 (t, r) + (2 + fc 2 ) U31 (t, r) = /„ (t, r) , 

9*tt33 (*, t) - 9fc 2 v 3 3 (t, t) = g u (t, t) , 

d t v 33 (t, t) + (2 + 9k 2 ) u 33 (t, t) = g v (t, r) . 



(36) 
(37) 



(38) 
(39) 
(40) 
(41) 



(42) 
(43) 



(44) 
(45) 
(46) 
(47) 



The last two equations do not give us relevant information for our purposes. On the contrary Eqs. (|4j) and (|4q) 
determine the evolution equation for the complex amplitude of the roll r. This comes from the fact that these equations 
contain resonant terms which yield divergent solutions unless a solvability condition is imposed. This is clearly seen 
by writing Eqs. (E3) and (Eq) in vector form: 



"3i (t,r) 




_ V31 (t, t) _ 







(2 + fc 2 ) 



fc 2 






"31 (*) t) 




' fu(t,T)~ 




_ v 31 (t, t) _ 


+ 


Jv(t,T)_ 



(48) 
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Upon diagonalizing this equation we obtain 

d tX (t,T)=-in X (t,T)+f(t,T), (49) 

(and its complex conjugate) where 

X (t, r) = (2 + k 2 ) u 31 (t, t) + iflv 31 (t, r) , (50) 

/ (t, r) = (2 + k 2 ) j u (t, t) + iSlf v (t, t) . (51) 

Eq. ( (49|) can only be solved if the driving term / does not contain elements oscillating as exp (— iClt). Upon substituting 
Eqs. (|38| ) and (39) into Eq. (|5l|), and making use of Eqs. (|22]), (p3|) and (32)-j35|), the solvability condition is found 
to be 

dr (l , , a 2 k2 2iuj 2 T- .3 + 5fc 2 2 

— = -72 (1 + fc jr + z— — e 2 r-z - \r\ r. (52) 

Finally we define a new roll complex amplitude 

R(t) = er (t) e-' l ^- T , (53) 
and turn back to the original parameters by undoing the scalings (|l^) with the result 

f = - [ 7 (1 + k 2 ) + i (u SI)] R + l ^R z^±^ |i?| 2 R. (54) 

Eq. ( |5~4| ) is the searched roll amplitude equation. It is a Landau equation with broken phase symmetry (note the 
presence of the linear term proportional to R). The roll solution ( |l5| ) can be written in terms of R making use of Eqs. 
@, @, Q and @. To the leading order the roll reads 

il> (x, t) = e - t(a/ul) sin(2wt) [1 + w (t) cos (kx)} , (55) 
w (t) = (1 - n/k 2 ) R (t) e lut + (1 + n/k 2 ) R (t) e- luJt . (56) 

Let us finally note that a straightforward linear stability analysis of the trivial solution R = (hence w = 0) of Eq. 
( |54| ) yields the same neutral stability curve given in (p^). 
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